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Abstract-  Output  signal-to-noise  ratios  for  a  multichannel 
decision  feedback  equalizer  were  measured  from 
experimental  BPSK  data  as  a  function  of  the  receiver  spacing, 
number  of  receivers  used,  and  array  aperture  and  reported  in 
a  paper  in  OCEANS  2005  by  the  author.  This  paper 
discusses  the  theory  and  reports  numerical  modeling  results 
to  compare  with  the  data. 

L  INTRODUCTION 

Noise  and  inter-symbol  interference  (ISI)  are  two 
primary  causes  of  acoustic  communication  bit  errors. 
Spatially  separated  multiple  receivers  can  be  used  to  (1) 
enhance  the  signal-to-noise  ratio  (SNR)  as  in  array 
beamforming  (or  equivalently  suppression  of  the  noise 
relative  to  the  signal)  and  (2)  suppress  the  ISI  (or 
equivalently  enhance  the  signal-to-interference  ratio,  SIR) 
using  diversity  combining  algorithms.  The  spatial 
processing  gain  can  be  illustrated  using  the  passive-phase- 
conjugation  (or  passive  time-reversal)  algorithm.  The 
received  data  can  be  expressed  in  terms  of  the  transmitted 
symbols  as  follows: 

L 

Vk  =  ^  K^k-n+rlk  (1) 

n= 0 

where  k  is  sample  index,  {/?„}  is  a  discrete  representation 
of  the  band  limited  impulse  response  timet  ion  of  length 
L+ 1,  and  {v*},  {%}  are  discrete  received  signal  and  noise 
sequences.  Convolving  Eq.  (1)  with  the  complex 
conjugate  of  the  channel  impulse  response  function  at  time 
t0,  and  summing  over  the  receiver  channels,  with  index  j, 
one  obtains  after  some  manipulations  in  summation  indices 
[1], 

L 

yk  =  Z  h*vk-t 

,  (2) 

L  L  L  v  ' 

—  X(J k  Xn^k-n  X  J k+n  K  *lk-l 

n= 1  n= 1  /=0 


M  L-n 

where  xn  =  is  the  discrete 

7=1  /=i 

representation  of  the  auto-correlation  of  the  channel 
impulse  response  function,  referred  to  as  the  Q  function 
[2],  (*  denotes  complex  conjugation.)  The  first  term  on 
the  right  hand  side  of  Eq.  (2)  shows  that  the  signal  is 
enhanced  using  multiple  channels.  For  a  time-invariant 

communication  channel  A.  ,{t)  =  A.  ;(/0)  ,  the  (matched 
field)  signal  gain  for  x0  is  ~101ogM  for  M  channels 
assuming  hj  —  .  The  second  and  third  term  show  the 

origin  of  ISI  due  to  non-zero  Xn  for  n  =  1,..  L.  The  x„  for  n 

V  0  are  the  sidelobes  of  the  Q  function  [2].  Previous  work 
has  shown  that  the  ISI  is  reduced  using  multiple  channels 
[2],  The  reason  is  that  the  sidelobes  appear  at  different 
positions  and  are  reduced  by  approximately  1/M  when 
averaged  over  M  channels. 

For  a  time-varying  channel,  h-  (t)  hj  (t0)  .  The 

matched-field  signal  gain  suffers  a  loss  as  a  function  of 
time.  Figure  1  shows  the  measured  channel  impulse 
responses  functions  for  intra-packet  correlations  and  the 
decrease  of  signal  gain  as  a  function  of  time  for  x o-  From 
the  physics  point  of  view,  signal  gain  depends  on  the 
signal  coherence.  The  loss  of  signal  temporal  coherence  is 
responsible  for  the  decrease  of  the  matched  field  signal 
gain  as  a  function  of  time.  From  the  signal  processing 
point  of  view,  the  signal  gain  degradation  is  caused  by  the 
tracking  error,  namely,  that  the  fixed  tap  coefficients 

hj  (/() )  fail  to  track  the  temporal  variation  of  the  channel. 

Note  that,  on  the  other  hand,  ISI  suppression  depend  on  the 
incoherent  or  independent  components  of  the  signal.  In 

fact,  X  for  n  ^  0  becomes  smaller  as  the  signal  becomes 

more  independent  (or  less  coherent)  between  the  receivers. 

This  paper  models  the  diversity  gain  arising  from  the  use 
of  diversity  combining  algorithms  applied  to  multiple 
receivers.  The  calculation  is  more  complicated  than  the 
passive  phase  conjugated  cased  discussed  above.  Classical 
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diversity  is  defined  for  receivers  receiving  independent 
information.  In  reality,  the  signals  on  spatially  separated 
receivers  are  never  totally  independent.  The 
corresponding  diversity  (for  partially  correlated  signals),  as 
measured  by  the  ability  of  the  receivers  to  suppress  ISI  is 
estimated  from  the  data.  The  equivalent  signal  gain  for  the 
multichannel  decision  feedback  equalizer  is  deduced  and 
compared  with  the  array  gain  of  a  conventional 
beamformer  applied  to  the  same  data. 


1  /N0.  One  finds  a  spatial  processing  gain  of  lOlogM  The 
model  assumes  a  frozen  ocean  and  an  ideal  solution  for  the 
equalizer  tap  coefficients.  Hence  one  obtains  a  theoretical 
processing  gain  equal  to  the  theoretical  array  gain  for 
conventional  beamforming.  In  the  real  world,  as  remarked 
above,  the  signal,  particularly  its  phase,  fluctuates  rapidly 
with  time,  resulting  in  limited  temporal  and  spatial 
coherence.  Hence  the  signal-to-noise  1/  N0  yield  less  than 
theoretical  lOlogM  gain. 

In  the  presence  of  ISI,  one  can  express  Eq.  (5)  as 


II.  OUTPUT  SNR  FOR  A  MULTI-CHANNEL 
DECISION  FEEDBACK  EQUALIZER 


The  output  (symbol)  SNR  for  communications  is 
defined  as  the  signal  to  noise-plus-interference  ratio 


Y  = - — - .  (3) 

K,  +isimt 

It  is  related  to  the  least  mean-square-error  (MSE)  by 

y=  1  y7’™ ,  (4) 


where  Jmin.  for  a  decision  feedback  equalizer  (DFE),  is 
given  below  [3,4],  assuming  that  the  channel  is  known: 


•Anin  -  N out  (1+ Iout  )> 


where 


f  rr t  7t/T  -*j  | 


Nn 


+  Nn 


(7) 

(8) 


for  additive  Gaussian  white  noise,  and 


=  exp 


1  +  ^Vq 
Q(0))  +  N{) 


\dco 


1.  (9) 


(2 n  7t  IT 

~z  f  ln 

-n ,T 


N n 


Q(co)  +  N0 


dco\,  (5) 


where  Q  is  the  normalized  spectrum  Q{cd)  —  Q{CO)  /  q0 , 

M 

where  Q{CO)  =  'ZQ,(a),  Qj {Cd)  -  H) (CO)Hj ( 0) ) , 

j= i 

with  j  being  the  channel  index,  H(a>)  is  the  Fourier 
transform  of  h(t)),  and 


ttIT 


do 


2k 


-71  IT 


J  Q(co)dco-^—  J  Qj  (co)d co  =  Mq0 


n!T  M 


-nIT  )= 1 

.2 


a -  2 

No  is  the  normalized  noise  Nn  =  — —  ,  where  (7~  is  the 


input  noise  level  (assuming  the  same  noise  level  for  all 
receivers).  We  note  that  qf)  is  the  mean  input  signal  level. 


and  N0  is  the  inverse  of  average  input  SNR  divided  byM. 
q0 ,  N0  and  Q((0)  are  dependent  on  the  number  of  channels. 
The  corresponding  values  for  M  channels  will  be  denoted 
as  q0M  ,  N0:m,  etc.  In  the  absence  of  ISI  ,  Q{CO )  =  1  one 
finds  Jmin  =  No/(l+No).  The  output  symbol  SNR  is  y  = 


Note  s{CO)  =  1  —  Q{CO )  ^  0  is  the  cause  of  ISI. 
We  note  for  <  1,  Eq.  (4)  can  be  written  as 


1 

NoJi+L„) 


(10) 


Comparing  Eq.  (10)  with  Eq.  (3),  one  finds  ISIou,=Nou,  Iout 
Equation  (10)  can  be  expressed  as 

OSNR(dB)  =  mog(r) 

=  -101og(M„„)-101og(l  +  /„„)’  (11) 

where  Nou,  is  in  essence  the  inverse  of  output  symbol 
energy  over  the  noise  and  1  +Ioui  represents  the  degradation 
due  to  ISI. 

Equation  (11)  can  be  stated  as  follows:  the  output 
(symbol)  SNR  can  be  improved  by  suppressing  the  noise 
relative  to  the  symbol  energy  (Nout)  and  the  interference 
relative  to  the  symbol  energy  (/„„,).  For  the  interference, 
note  that  =  0,  when  Q(co)  =  1.  In  other  words,  ISI 
comes  from  spectral  Q((0)  often  has  due  to  multipath 
interference  or  signal  fading.  The  spectral  nulls  in  Q( co)  are 
minimized  by  averaging  over  independent  channels  (which 
possess  nulls  at  different  frequencies).  This  is  the  reason 
that  spatial  diversity,  originally  conceived  to  combat 
temporal  fading,  can  be  used  to  mitigate  ISI.  In  a  random 
ocean,  the  acoustic  signals  encounter  (repeated)  scattering 
from  the  ocean  inhomogeneities  and/or  time-varying 
surfaces  (of  scales  smaller  than  the  spatial  coherence 
length).  The  spatially  uncorrelated  (independent) 
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scattering  components  of  the  acoustic  signals  help  to 
minimize  ISI. 

Equation  (11)  can  be  translated  into  an  expression  for 
the  spatial  processing  gain  defined  as  the  output  SNR  for 
M  channels  over  that  of  a  single  channel.  One  finds 


SPGm  {dB)  =  NSGm  +  ISGm  ,  (15) 

where 


NSGm  =10  log 


Nout(M) 

-  lOlogM  +  101og 


f  -  \ 

Clo,M 


V  ) 


+ofy 

do 


and 


(16) 


BG“=10lo8(i1S)’  (17) 

where  Nout{M)  and  Iou,{M)  are  given  in  Eqs.  (8)  and  (9) 
using  N0  and  Q  for  M  channels.  The  second  expression  in 
Eq.  (16)  assumes  a  high  input  SNR.  Equation  (15)  states 
that  spatial  processing  gain  (SPG)  is  the  sum  of  noise 
suppression  gain  (NSG)  and  ISI  suppression  gain  (ISG) 
expressed  in  dB. 


III.  EXPERIMENTAL  DATA 

In  a  previous  paper  [5],  we  reported  the  spatial 
processing  gain  as  a  function  of  number  of  receivers  using 
(1)  conventional  beamforming  and  (2)  a  multichannel  DFE 
jointly  with  a  phase-locked  loop  [1,2].  The  data  were  linear 
frequency  modulated  (LFM)  and  binary-phase-shift-keying 
(BPSK)  signals  from  the  ASCOTOl  experiment  which 
took  place  off  the  coast  of  New  England  [6].  Both  signals 
were  centered  at  1200  Hz  with  a  400  Hz  bandwidth.  The 
receiver  array  consisted  of  33  elements  uniformly  spaced 
at  d  =  0.5  m  covering  a  depth  span  of  50-66  m.  The  source 
was  deployed  4  m  above  the  bottom.  The  water  depth  was 
~100m.  The  source-receiver  range  was  approximately  10 
km.  The  received  signals  had  a  high  input  SNR  (>  15  dB) 
[5], 

We  processed  the  data  using  subarrays  of  different 
number  of  receivers,  with  different  receiver  separations:  D 
=  nd,  n  =  1 ...  10,  where  d  =  0.5  m.  Signal  gain  was 
measured  by  beamforming  the  various  sub-arrays  using  the 
LFM  signal.  To  obtain  array  gain,  we  assume,  for 
simplicity,  uncorrelated  noise  and  hence  a  noise  gain  of 
lOlogM,  where  M  is  the  number  of  receivers  of  the 
subarray.  The  AG  data  (see  [1,2])  can  be  fitted  with  AG  = 
A  log(l+L/D),  where  L  is  the  array  aperture,  and  A  is  a 
coefficient  which  decreases  linearly  with  increasing 
element  spacing:  A  lies  between  8  and  6  for  D  =  Id  to  lOd. 

The  output  symbol  SNR  for  the  multichannel  DFE  [7] 
was  measured  using  the  BPSK  signals  for  the  subarrays 
with  different  spacing  [5].  The  results  were  plotted  as  a 
function  of  the  array  aperture  as  shown  in  Fig.  2.  The  data 


can  be  fitted  by  the  following  formula:  output  SNR  = 
101og(L/p)  +  G,  where  G  =  8  for  D  =  Id,  2d,  3d,  4d;  G  = 
5.5  for  D  =  5d,  6d;  and  G  =  2.5  for  D  =  7d,  8d,  9d,10d. 
The  data  were  also  processed  as  a  function  of  receiver 
spacing  for  a  fixed  number  of  receivers  (M  =  4)  as  shown 
in  Fig.  3. 

The  above  result  suggests  that  for  a  given  array 
aperture,  optimal  performance  is  obtained  by  spacing  the 
element  at  D  <  p  (Fig.  2),  where  p  is  the  spatial  de- 
correlation  (or  coherence)  length  of  the  signal  p  =  4d.  To 
keep  the  number  of  elements  (and  computational 
complexity)  to  a  minimum,  one  should  set  the  element 
spacing  at  D  =  p.  Note  that  increasing  the  number  of 
receivers  using  smaller  separations  does  not  necessarily 
improve  the  output  SNR  (Fig.  2).  The  result  for  a  fixed 
numbet  of  receivers  (Fig.  3)  also  suggests  that  the  optimal 
receiver  spacing  is  approximately  the  signal  de-correlation 
length. 
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Fig.  1  Channel  Impulse  response  as  a  function  of  delay 
time  and  geotime  (upper  figure).  Signal  gain  degradation 
as  a  function  of  geotime  (lower  figure). 

To  compare  and  interpret  the  experimental  results,  we 
calculate  numerically  the  SPG  using  Eqs.  (15-17).  The 
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input  is  the  channel  transfer  function,  which  was  estimated 
using  the  LFM  signals.  We  note  that  Eqs.  (15-17)  assume 
a  time-invariant  channel  which  is  known  to  the  receiver. 
In  other  words,  the  tap  coefficients  are  well-matched  to  the 
channel.  In  the  real  world,  there  is  a  mismatch  between 
the  tap  coefficients  and  the  channel,  and  hence  the 
measured  SPG  will  be  less  than  the  theoretical  SPG. 

The  SPG,  as  given  in  Eq.  (15)  is  the  sum  of  noise 
suppression  gain  (NSG)  and  ISI  suppression  gain  (ISG). 
The  ISG  is  modeled  using  Eq.  (17).  The  results  for  the 
various  subarrays  are  shown  in  Fig.  4.  One  finds  that  the 
ISG  is  practically  only  a  function  of  the  array  aperture  and 
is  nearly  independent  of  the  element  spacing.  For  the 
various  array  configurations  mentioned  above  with 
different  element  spacing,  the  ISG  can  be  fitted  with  the 
following  formula:  ISG  =  101og(l+L/p)  -  1.  This  result 
suggests  that  populating  the  array  with  more  elements  does 
not  always  improve  the  ISG  because  the  additional 
receivers  are  NOT  receiving  new  (independent) 
information.  This  is  not  necessarily  true  for  the  SPG  (or 
the  output  SNR). 


Fig.  2  Output  SNR  as  a  function  of  the  array  aperture 
with  different  receiver  separations.  Thick  black  curves  are 
empirical  fit  to  the  data  (see  text). 

One  finds  that  the  NSG  based  on  the  theory  is  nearly 
equal  to  the  theoretical  maximum  of  lOlogM  = 
101og(l+L/D).  This  is  expected  based  on  Eq.  (16). 
Combining  ISG  with  NSG,  one  obtains  a  theoretical  SPD 
as:  SPDth  =  101og(l+L/D)  +101og(l+L/p)  -  1.  In  contrast, 
the  output  SNR  data  in  Sec.  II  produces:  SPDData  = 
101og(L/p)  +  G  -  OSNR(l),  where  the  last  term  is  the 
output  SNR  for  one  channel. 

The  theoretical  SPD  result  suggests  that  for  a  fixed 
array  aperture  L,  optimal  SPD  is  obtained  by  spacing 
element  at  half  wavelength,  D  =  d  since  SPDlh  increases 
with  decreasing  D  (NSG  is  maximum  at  D  =  d).  In 
contrast,  as  noted  above,  the  data  suggest  that  SPD  for  D  = 
Id,  2d,  3d,  and  4d  are  approximately  the  same.  The  reason 
is  that  the  NSG  in  the  real  world  is  different  and  smaller 


than  the  theoretical  value,  lOlogM.  Note  that  in  the  real 
world,  the  signals  are  time-varying.  Note  further  that  the 
enhancement  of  the  signal  over  noise  (the  NSG)  requires 
coherent  processing  of  the  signal.  The  tap  coefficients, 
estimated  based  on  the  least  mean  square  error  criterion, 
could  not  match  the  temporal  variation  of  the  channel. 
Consequently  the  NSG  of  the  equalizer  is  less  than  the 
theoretical  value. 

The  NSG  in  a  time-varying  channel  is  difficult  to 
model  analytically.  Given  the  measured  output  SNR  from 
data,  one  can  estimate  the  NSG  by  subtracting  the  ISG 
from  SPDData.  Note  that  ISI  involves  the  independent 
(incoherent)  component  of  the  signal  and  can  be  reliably 
estimated  from  the  spectrum  of  the  signal.  Using  data  from 
Fig.  2  and  4,  we  obtain:  NSGdata  =  -101og(l+p/L)  +  Q,  Q  = 
G  +  C  -  OSNR(l). 

Using  the  formula  for  ISG  and  NSGdata,  one  can 
calculate  the  SPD  for  an  array  of  four  receivers  and  then 
the  output  SNR  for  four  channels.  The  modeled  output 
SNR  for  four  receivers  is  shown  in  Fig.  3  as  the  dotted 
line.  It  is  in  good  agreement  with  the  data  reported  above. 

For  conventional  beamforming  the  measure  of 
performance  is  the  array  gain.  The  measured  array  gain  for 
the  ASCOTOl  data  was  given  above.  For  the  multichannel 
combining  algorithm,  the  measure  of  performance  is 
diversity  gain.  What  is  diversity  gain?  There  are  two 
definitions.  Diversity  gain  has  been  used  loosely  in  the 
literature  to  refer  to  the  improvement  in  the  symbol  SNR 
through  the  use  of  a  diversity-combining  algorithm  applied 
to  widely  spaced  receivers.  This  was  the  definition  used  in 
[5]  and  is  referred  to  above  as  the  spatial  processing  gain. 
Classical  diversity  requires  that  the  signals  on  the  different 
receivers  fade  independently  and  are  thus  uncorrelated 
between  the  receivers.  For  this  case,  the  signal  is  like  noise 
and  cannot  be  enhanced  over  noise,  or  equivalently  no 
NSG  is  obtained.  The  performance  improvement  comes 
primarily  from  the  suppression  of  ISI.  Thus  classical 
diversity  gain  is  more  in  line  with  the  ISG.  This  is  the 
second  definition  of  diversity  gain.  Note,  however  that 
classical  diversity  does  not  exist  in  the  real  world,  since  the 
signals  are  never  totally  uncorrelated.  For  partially 
correlated  signal,  the  ISG  is  the  “available”  diversity  gain. 
The  formula  for  it  is  given  above. 

IV.  SUMMARY 

In  this  paper,  we  modeled  the  spatial  processing  gain 
for  a  multichannel  DFE  in  terms  of  the  signal  to  noise 
enhancement  (NSG)  and  signal  to  interference 
enhancement  (ISG)  and  compare  the  results  with  the 
measurements  reported  in  a  previous  paper.  The  available 
diversity  gain  for  partially  correlated  signals  is  derived. 
The  spatial  processing  gain  for  a  multichannel  DFE  are 
fitted  in  terms  of  the  signal  aperture  and  signal  coherence 
length  and  compared  with  the  array  gain  of  conventional 
beamforming.  Details  can  be  found  in  [8]. 
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Fig.  3  Output  SNR  as  a  function  of  receiver  spacing  for 
an  array  of  four  receivers.  Solid  curve:  data.  Dashed 
curve:  model. 
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Fig.  4  Modeled  ISG  as  a  function  of  the  array  aperture 
with  different  receiver  separations.  Thick  curve  is  an 
empirical  fit  to  the  ISG  as  discussed  in  text.  Same  legend 
as  in  Fig.  2. 


5 


